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Abstract 



K^ i We introduce a new formulation of the diquark breakins; mechanism which 

(N 

en ■ describes with no free parameters the huge nuclear stopping observed in cen- 

(N 

psj ' tral nucleus-nucleus collisions. Supplemented, in the dual parton model, with 

o : 

Q"^ . strings originating from diquark-antidiquark pairs from the nucleon sea, it 

gives a substantial increase of hyperon and antihyperon yields from pPb to 

Qh! central PbPb collisions. Compared to data, this increase is slightly underes- 

Mh. timated for A's and H's but is five times too small for O's. Introducing final 

"^ ' state interactions, a reasonable description of all baryon and antibaryon yields 

is achieved. 
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Enhancement of strange baryon and anti-baryon production in central nucleus- 
nucleus collisions has been observed experimentally [Up. In particular, a spec- 
tacular enhancement of H's and fi's from pPb to central PbPb collisions has been 
measured |^ . This phenomenon has been proposed as a signal of quark gluon plasma 
formation |^. On the other hand explanations based on string models have also been 
proposed |^. In ref. enhancement of A's and A's was described in the dual parton 
model (DPM) with final state interaction. An important drawback of [^ is the 
lack of nucleon stopping -which would increase the net baryon yield at mid rapidi- 
ties. As a consequence a very strong effect due to final state interaction was needed, 
as well as a rather large value of the fraction of strange quarks in the nucleon sea. A 
mechanism of nuclear stopping based on diquark breaking (DB) has been known for 
a long time 0-0. Recently, it has been introduced in nucleus-nucleus interactions 
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The purpose of this paper is threefold. First, we propose a new formulation 
of the DB stopping mechanism of |T^-[|ri|], which allows to describe net baryon 
stopping without any dynamical parameter governing the DB probability. Second, 
we show that the DPM with DB and strings with diquarks-antidiquark pairs (d-d) 
from the sea at their ends, allows to describe most of the observed yields of A, A, S 
and S -but significantly underestimates the production of fi's and O's. Third, using 
the corresponding baryon densities as initial conditions for final state interaction we 
achieve a reasonable description of all baryon and antibaryon production^. 

Baryon stopping . We describe baryon stopping via the mechanism shown in Fig. 1 



-the same as Fig. 5 of ||rT|. However, contrary to [|Tl|] we assume here that there is 



no dynamical suppression associated to DB -i.e. the weight of the DB diagram of 



^In this paper we concentrate ourselves on the rapidity density. For a calculation of the trans- 
verse mass spectra in DPM see [n2| . 



Fig. 1 is the same as the conventional DP one of Fig. 2. More precisely, in the case of 
a single inelastic collision per nucleon (two string configuration), two valence quarks 
of the nucleon are sitting together at a string end. We assume that they form a 
sort of compact object which hadronizes in the conventional way (DP mechanism). 
However, in the presence of several inelastic collisions we assume that the net baryon 
number (string junction) has equal probability of finding itself in any of them. If it 
follows the valence diquark (Fig. 2) we have the conventional DP situation. In the 
other cases it joins with sea quarks producing a much slower baryon (DB mechanism 
of Fig. 1). The rapidity distribution of the net baryon production AB = B — B in 
AA collisions is then given by 



iy) 






(1) 



dy n 

Here ua is the average number of participants in each nucleus and n the average 
number of collisions. dNDp/dy is given by the conventional DP hadronization me- 
chanism and dNoB/dy is the rapidity density given by the DB component -which 
behaves [0, ^, |IU| as exp[— |||/ab— yAfaa;|]-| Both distributions depend, due to energy- 



conservation, on the average number of collisions per nucleon ti/ha ~ n/fiA- Note 
that for n = fiAi only the DP component is present and we recover the usual scaling 
in the number of participants. Baryon number conservation implies that both DP 
and DB rapidity distributions are normalized in such a way that their integral over 



y is equal to two. Following [|ri[] we have 



■''This is the behaviour in the original work 0] and is used by aU authors at present energies 
llOl , [ll[ . However, a fiat behaviour in y at high energies has been proposed in [g|, Q. This has 
very important consequences for the net baryon number aX y* ^ Q at LHC |0]. A new component 
consisting on a diquark which contains a quark from the nucleon sea also has been introduced in 
|13| . However, in this case the diquark hadronizes as a valence diquark producing baryons mainly 
in the fragmentation regions. 



-^(y) = Cnl^A {ZT{^ - Z+V-'^- + Z'J\l - Z.Y'-/-a) (2) 

where Z± = exp{±y — i/max)- The same shape will be used for all baryon species. 



Two different ansatzs for the power 7 were proposed in |TT[] -which give for central 
PbPb collisions (n/fiA ~ 4) a value of dN^^/dy at y* ~ in the range 0.5 to 0.60. 
In order to reproduce the data of the NA49 collaboration [|l| for B-B at y* = 0, we 
need a value 0.48. Our analysis will concentrate on mid-rapidities and this value 
will be used throughout this paper. 

Eq. ([1|) also applies to pp collisions putting n = ua = I. However, this ap- 
plies only at low energies where the two string configuration dominates. At high 
energies the second term in (|I|) is present and increases the net baryon stopping. 
Unfortunately, the only available data at ISR do not allow to either confirm or rule 
out the presence of the DB component. However, the effect of this component is 
quite important in pA collisions where the pPb data of ^ confirm its presence (see 
below). In DPM, the extension of eq. (|l]) to pA collisions is straightforward. A 
good approximation to the exact formula, valid at y* = 0, is 

Here n = A app/apA and {dN^^ /dy)n is given by eq. (P with ua = 1. Eq. (H) 
is self-evident. On the nucleus side we have n diquarks which hadronizes as in pp 
collisions (the average number of collisions per nucleon at present energies is close 
to one in both cases). On the proton side, we have one diquark which hadronizes 
as in AA collisions. (The average number of collisions of this nucleon in pPb and in 
central PbPb collisions is practically the same). Eq. (|^) is very useful since it gives 
the rapidity density of the net baryon in pPb, in terms of the same quantities which 
appear in PbPb (eq. (|l])) plus the corresponding one for pp. 



B-B pair production . Following p, |T^ the rapidity distributions of antibaryons 
in AA collisions is given by 



Here dNstring/dy denotes the conventional pair production resulting from d-d pro- 
duction in the string breaking process. At present energies the q-q strings have too 
little energy for baryon pair production ; practically all their production occurs in 
qq-q strings, hence the scaling in the number of participants. This gives a much too 
small antibaryon yield. For this reason we have introduced p, |12[ the second term 



in (^ which corresponds to pair production in strings with a,d oi d from the nucleon 
sea at one of their ends. The average number of d-d pairs scales as n-fiA- Note that 
for dNstring/dy ~ dNsea/dy, the scaling in the number of participants is converted 
into a scaling in the number of collisions. 

For pA collisions the equivalent of expression (^, vahd at j/* = 0, is 

Here n = A app/apA and {dN^ /dy)n is given by eq. (H) with ua = 1. 

Numerical results . Let us concentrate first on the non-conventional contributions 
(second terms in eqs. (^^ and (@-^) -which are numerically very important since 
they contain the factor n-fiA- Actually, dNf^^/dy has been computed in [p. However, 
as stated there, there is a huge uncertainty on its absolute normalization. In view 
of that, we are going to determine it from pA data. Since data for multi-strange 
hyperons are only available in the rapidity range —0.5 < y* < 0.5 we shall restrict 
our analysis to this region. 



Turning to the relative yields of the different baryon species, we note that in 
dNf^^/dy the baryons and antibaryons are produced from three sea quarks. There- 
fore the factors controlling the relative yields are : 4L^, 4L^, 12L^5', SLS*^, "iLS^ and 
S^ for p, n, A + S, S°, S^ and fi, respectively. Moreover, we use for both baryons 
and antibaryons the conventional relation S+ + S^ = 0.6A, which is obtained after 
resonance decay. We take L = 0.435 and S = 0.13. In this way the strangeness 
suppression factor has the conventional value S/L = 0.3 and the sum of all baryon 
yields add up to (2L + S)'^ = 1. However, this fixes only the relative normalization 
of the various (anti) baryon species -and there is one free parameter which provides 
the absolute normalization which is fixed by the data. The obtained values are given 
in Table I. 

The situation is quite different for the new component dN^^/dy. Here net 
baryon conservation fixes the absolute normalization and therefore there is in prin- 
ciple no free parameter -since the strangeness suppression factors are the same as 
above. However, it turns out that proceeding in this way the strange baryon yields 
in pA are overestimated -especially S's and fi's. A possible explanation of this fact 
is that, at present energies, it may happen that the baryon is not formed out of three 
sea quarks as in Fig. 1, but, due to lack of phase space, the valence quark at the 
string end is picked up instead^. In this case, the relative yields of A's of S's will be 
reduced and that of fi's will be zero. In any case, we take here a phenomenological 
attitude. We determine the strange net baryon yields at y* = from the data. We 
have thus three free parameters A-A, S^-H+ and Q-Cl -plus a forth parameter for 
the absolute normalization of antibaryons as discussed above. We arrive in this way 
to the values listed in Table I. 



''Nevertheless, the factors given above should apply at very high energy. They give a ratio 
S/A + S -- 0.3 in agreement with Tevatron (wl and SPS coUider 111 data. 



Finally, the conventional components dNcp/dy and dNf^^^^g/dy have been ob- 
tained in two different ways : using the formalism of ref. P] and using the QGSM 
T6| , p^ Monte Carlo without string fusion]^. In the former case, the obtained values 



are given in Table I. It should be noted that the contribution of these conventional 
components to the final results is significantly smaller than the one of the new 
components, since the latter are proportional to n-fiA- 

With the values given in Tables I we obtain the results shown in Figs. 3 and 4. 
The pA data are well reproduced for all baryon species. In PbPb, the p + p and 
most of the A + A yield are reproduced. However, the S^ + S+ yield is somewhat 
underestimated and the fi + one is five times smaller than the measured on^. 

In an attempt to solve this problem, we use our expressions (|l]) and (^ as initial 
baryon densities for final state interaction. We consider the following channels 

nN ^ KA , ttN ^ KE , nA ^ KE , vrS ^ iTS , nE ^ Kil , (6) 

plus the corresponding ones for antiparticle^^ The final state interaction is treated 

^The parameters of this Monte Carlo, which produced too many strange baryons in pp coUisions 
Jl7], have been tuned up, and its results in pp collisions are now in agreement with the calculations 

of ^ and with the experimental data (see Table 1 in [g|). 

"a qualitatively similar result has been obtained in ||l8|] . This paper, which was appeared after 
completion of the present work, is based on a different realization of baryon stopping and uses a 
junction-antijunction exchange mechanism for baryon pair production. It does not include final 

state interaction. 

**As a technical point, one should notice that some of the charge combinations in (jm are not 

possible or have negligeably small cross-sections. The latter correspond to quark diagrams with 

three quark lines (baryon exchange) in the i-channel. They have not been included. All reactions 

corresponding to quark diagrams with light quark annihilation and s-s production |j have been 

included with identical cross-section. We have neglected all strangeness exchange reactions ttA -^ 

KN, etc. Although, at threshold, the cross-sections are larger than those of the non-strangeness 

exchange reactions (^, this is no longer the case for the averaged cross-sections < a > (see Q). 

Channels (m are thus dominant due to the relations pjv > PA > Ps > Pn between baryon densities. 



in the same way as in ref. [|1^, ^. The products of densities in (H) are taken at 
fixed values of b and s. For the baryon densities, the dependence on these variables 
is contained in the geometrical factors ua and n in eqs. (|I|) and (^). They are 
computed in the Glauber model with Saxon- Woods profiles. The pion densities and 



the initial and final times are the same as in [|T^, ^. The only difference resides in 
the fact that, in the case of J/ip suppression, there is a single channel and the J/ip 
appears both in the initial and final state. In this case, the differential equations 
for the gain and loss of particles can be solved analytically and a simple exponential 
form is obtained. In the present case, the gain and loss differential equations for 
the solution of the coupled channels (|) has to be obtained numerically. We have 
only one free parameter : the cross-section < cr > averaged over the momentum 
distribution of the incoming particles which is taken to be the same for all reactions 
(the relative velocity between the interacting particles is included in < cr >). We 
take < cr >= 0.14 fm. (A comparable value has been obtained in ||^). It is clear 
that the relative enhancement of the various baryon and antibaryon yields resulting 
from the interactions (j^) will be the largest for f2's and the smallest for A's. This 
is due to the fact that A's are produced in some reactions and destroyed in others. 
On the contrary fi's are produced and never destroyed. The results are shown in 
Figs. 3 and 4. A reasonable description of all baryon and antibaryon is obtained. 
However, the ratio Q + Q/E^ + S+ in PbPb collisions is slightly underestimated. 

In conclusion, we have introduced a new formulation of the diquark breaking 
mechanism which describes baryon stopping with a diquark breaking probability 
equal to one. Incorporated in DPM, together with strings originating from diquark- 
antidiquark pairs from the nucleon sea and final state interactions we have obtained 
a reasonable description of all baryon and antibaryon yields at mid-rapidities in pp, 



pPb and PbPb collisions. 
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Figure Captions 

Fig. 1 : Diquark breaking component for nucleon-nucleus scattering and two 
inelastic collisions. 

Fig. 2 : DP component for nucleon-nucleus scattering for two inelastic collisions. 

Fig. 3 : Baryon yields for pPb and PbPb collisions at 160 GeV per nucleon 
calculated from eqs. (^^ using the values in Table 1 before (dotted line) 
and after (full line) final state interaction. The dashed lines are the results 
with final state interaction using the Monte-Carlo results |1^ |17] for the 



conventional DP and string breaking components. The data are from the WA 
97 collaboration 0. 

Fig. 4 : Same as in Fig. 3 for the ratios B/B. 
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Table 





P 


A 


S 


n 


dN^E/dy 
dNpJdy 

dN^E/dy 
dNl^Jdy 

dN^^Vdy 
dN^i/dy 


1.87 X 10-1 
4.36 X 10-3 
6.90 X 10-2 
8.50 X 10-3 
3.60 X 10-2 
1.70 X 10-2 


6.22 X 10-2 
2.44 X 10-3 
1.40 X 10-2 
2.83 X 10-3 
7.00 X 10-3 
5.65 X 10-3 


3.26 X 10-3 
2.92 X 10-^ 


1.65 X lO-'^ 


3.30 X 10-^ 


2.38 X 10-^ 
2.91 X 10-^ 


5.05 X 10-6 


1.01 X 10-5 



Table I : Values of the rapidity densities at y* = in eqs. (|T],^|D. The values in 
the first lower rows are for central PbPb collisions (n^ = 178, n = 858) and those 
in the last two rows for pp. Three numbers in the first row and one in the second 
one have been adjusted to the data (see main text). 
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Figure 1 
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Figure 2 
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Figure 3 



c 

> 

0) 



(f) 

c 

o 



10 r- 



>^ 1 



^ 



-2 



^ 



-4 




A+A 



+ 



10 



N 



part ' 



16 



Figure 4 
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